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ABSTRACT 

We present extensive radio and X-ray observations of the nearby Type Ic SN 2007gr in NGC 1058 obtained 
with the Very Large Array and the Chandra X-ray Observatory and spanning 5 to 150 days after explosion. 
Through our detailed modeling of these data, we estimate the properties of the blastwave and the circumstellar 
environment. We find evidence for a freely-expanding and non-relativistic explosion with an average blast- 
wave velocity, v « 0.2c, and a total internal energy for the radio emitting material of £ « 2 x lO"** erg assuming 
equipartition of energy between electrons and magnetic fields (6^ = e^ = 0.1). The temporal and spectral evo- 
lution of the radio emission points to a stellar wind-blown environment shaped by a steady progenitor mass 
loss rate of M sa 6 x 10"^ M0 yr"' (wind velocity, Vh, = 10^ km s~'). These parameters are fully consistent 
with those inferred for other SNe Ibc and are in l ine with the e xpect ations for an ordinary, homologous SN 
explosion. Our results are at odds with those of Paragi et al. I (12010 ) who recently reported evidence for a 
relativistic blastwave in SN 2007 gr based on their claim that the radio emission was resolved away in a low 
signal-to-noise Very Long Baseline Interferometry (VLBl) observation. Here we show that the exotic physical 
scenarios required to explain the claimed relativistic velocity - extreme departures from equipartition and/or a 
highly collimated outflow - are excluded by our detailed Very Large Array radio observations. Moreover, we 
present an independent analysis of the VLBl data and propose that a modest loss of phase coherence provides 
a more natural explanation for the apparent flux density loss which is evident on both short and long baselines. 
We conclude that SN 2007 gr is an ordinary Type Ibc supernova. 

Subject headings: supernovae: specific (SN2007gr) 



1. INTRODUCTION 

Nearly 25 years have elapsed since Type Ibc supernovae 
(SNe Ibc) w ere first recognize d as a distinct class of cosmic 
explosions jElias et al. 1 119851; iFilip penko & Sargent 19851 
IWheeler & Levre ault 19851). They are now commonly under- 
stood to mark the gravitational core-collapse of massive stars 
from which an explosion launches ejecta to bulk velocities of 
- 10,000 km s"' (e.g., Filippenko' 1997). Motivated by the 
observed lack of hydrogen (and often helium) in their optical 
spectra, the favored progenitors of SNe Ibc are Wolf-Rayet 
stars that hav e ejecte d their massive envelope through strong 
stellar winds (iBegelm an & Sarazin 198 a) or inter action with 
a binary companion (Podsiadlows ki. Joss & Hsuiri992) . Re- 
cently, the observational realization that most long-duration 
gamma-ray bursts (LGRBs) are accompanied by SNe Ibc has 
fuel ed a new surge of interest in these massive star explosions 
(see lWooslev & Bloomll2006l and references therein). 

The LGRB-SN connection and the relative rates of these 
events implies that LGRBs are a rare subclass of SNe 
Ibc, distinguished by the production of a relativistic out- 
flow that decouples from the bulk SN ejecta during ex- 
plosion. For many LGRBs, the energy associated with 

' Harvard-Smithsonian Center for Astrophysics, 60 Gai'den St., Cam- 
bridge, MA 02138, USA 

- Max-PIanck-Institute for Radio Astronomy, Auf dem Hugel 69, 53121 
Bonn, Germany 

^ Raymond and Beverly Sackler School of Physics & Astronomy, Tel Aviv 
University, Tel Aviv 69978, Israel 

'^ University of Virginia, Astronomy Department, Charlottesville, VA 
22904, USA 

^ Hartebeesthoek Radio Observatory, P.O. Box 443, Krugersdorp, 1740, 
South Africa 

^ Dept. of Physics and Astronomy, York University, Toronto, M3J 1P3, 
Ontario, Canada 



^52 



the relativistic blastwave, E « 10 - 10 erg, is compara- 
ble to that of the non-relativistic SN ejecta, £sn ~ 10^' - 
10^^ erg (e.g GRB 030329/SN2003dh;jBerger ef a/ 
iHjorfli et g/. I I2003t iMazzaU et al. I I200I iFrail et alT 
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2005h . 



This duality challenges the standard SN mechanism in which 
a homologous explosion couples at most 0.01% of the total 
energ y to mildly-relativistic ejecta (iTan. Matzner & McKed 
l2001h . It is therefore generally accepted that a "cen- 
tral en gine" — an accreting and rapidly-ro tating compact 
object (MacFadven. Wooslev & Heger' '2 0011) or a magnetar 
(e.g., Thompson, Chang & Ouataert 20041) — is additionally 
required to power the energetic and relativistic blastwaves of 
LGRBs. The critical question is whether ordinary SNe Ibc 
also harbor central engines, and in turn, the ability to produce 
even weak relativistic outflows. 

Identifying engine-driven, relativistic explosions requires 
direct measurements of the blastwave velocity and energy. 
While optical SN emission predomi nantly traces the Nickel- 
56 synthesized in the bulk SN ejecta (lArnettlll982h . radio and 
X-ray observations directly probe the synchrotron radiation 
produced as the blast wave sho cks material in the circumstel- 
lar medium (CSM; Chevalienll982{) . Over the past decade, 
dedicated radio studies of SNe Ibc have consistently pointed 
to blastwave velocities of just (v) « 0.15c and associated 
energies, (E) w 10"*^ erg, in agreement with the expecta- 
tions for an ordinary core-collapse explosion (Chevali^ 
1998; Berger, Kulkarni & Chevalier 2002; Bergerefa, 
2003; Soderber g ef g/. ■ ^2005, 200 6; .Chevalier & Fransso 
2006; Soderbergfl^Ol [Soderberg efa/J |2008t). Recently, 
luminous radio emission was detected from the broad-lined 
Type Ic SN2009bb pointing to an energetic (E > 10'*'' erg) 
and rel ativistic (v > 0.85c) o utflow powered by a central 
engine (ISoderberg et al. l[2010l) . This discovery marked the 






first relativistic SN explosion identified without a detected 
gamma-ray counterpart. Together with the results of on-going 
radio supernova studies, this result independently confirmed 
that the fraction of SNe Ibc with engi ne-driven rela t ivistic 
outflows is exceedingly small, ^ 1 % (iBerger et al. 1 120031: 
ISoderberg et al. 2006. 2010) 

Against this backdrop, iParagi et oTI ( 1201 Ol) recently 
claimed evidence for a mildly-relativistic blastwave (v > 
0.6c) from the otherwise ordinary Type Ic SN2007gr based 
on a low signal-to-noise Very Long Baseline Interferome- 
try (VLBI) radio observation. Here we present an indepen- 
dent analysis of these VLBI data and report in detail on our 
extensive multi-frequency radio and X-ray observations of 
SN 2007gr obtained with the Very Large Array^ ( VLA) and 
the Chandra X-ray Observatory (C XO). 

In contrast to the hypothesis of iParagi et al. I ( 120101) . our 
multi-frequency radio data (j/obs = 1.4-22 GHz) and early X- 
ray observation reveal that the SN emission is best explained 
by an ordinary, non-relativistic blastwave with an average ve- 
locity, V w 0.2c, and an associated energy, £ w 2 x 10 erg. 
This result stems from the fact that we directly measure the 
properties of the radio spectral peak across multiple epochs. 
We further demonstrate that a conservative interpretation of 
the VLBI data is fully consistent with non-relativistic SN ex- 
pansion. Moreover, we argue that the exotic physical condi- 
tions required by a mildly relativistic outflow are ruled out by 
our VLA observations of the temporal and spectral evolution 
of the radio emission. Overall, we find that the SN 2007gr 
blastwave properties are analogous to those of other Type Ibc 
supernovae and starkly dissimilar from the relativistic jets of 
LGRBs. 

2. OBSERVATIONS 

SN2007gr was discovered by the Katzm ann Automatic 
Imaging Telescope on 2007 Aug 15.51 UT ([Ma dison & Li 
12007.) in NGC 1058 (d « 9.3 Mpc; ISilbermann7f a/. ..1996.) . 
An early spectrum of the supernova indicated a Type Ibc 
classification based on the preliminary identification of He I 
features (Chornockef a/. 2007) although it has been argued 
that the lines are better matche d by C II motivatin g the re- 
classification as a Type Ic SN (IValenti et al. 112008 1). Based 
on the non-dete ction of the SN in pr e-discovery images from 
Aug 10.44 UT Madison &Til (|2007|) . we adopt an explosion 
date of 2007 Aug 13 ± 2 days. 

2.1. Very Large Array 

We observed SN 2007gr with the VLA beginning on Aug 
17.41 UT (At « 4 days after explosion) as part of the Very 
Large Arra y Intensive Stud y of Naked Supernovae project 
(VISioNS; ISoderberd l2007h . We detected a radio source 
coincident with the optical position with a flux density of 
/V = 610 ±40 ^l]y atiy = 8.46 GHz (Soderberg 2007). We 
subsequently initiated an intense VLA follow-up campaign to 
study the temporal and spectral evolution of the radio emis- 
sion. 

Radio data were coflected at 1.43, 4.86, 8.46, and 22.5 GHz 
between 2007 August and 2008 January (Table [T]!. The bulk 
of our data were obtained in the highest resolution VLA ar- 
ray (A array), and all observations were taken in the standard 

' The Very Large An'ay and Very Long Baseline An'ay are operated by the 
National Radio Astronomy Observatory, a facility of the National Science 
Foundation operated under cooperative agreement by Associated Universi- 
ties, Inc. 



continuum observing mode with a bandwidth of 2 x 50 MHz. 
At 22.5 GHz we included reference pointing scans to correct 
for the systematic 10-20 asec pointing errors of the VLA an- 
tennas. We used 3C48 (J0137H-331) for flux calibration, while 
phase referencing was performed using calibrator J0254-I-395. 
Data were reduced using standard packages within the Astro- 
nomical Image Processing System (AIPS). We fit a Gaussian 
model to the radio SN in each observation to measure the inte- 
grated flux density, and for our 8.46 GHz data set (our highest 
resolution observations) we further report the associated cen- 
troid position (Table[Tli. 

The statistical positional errors we infer from each epoch 
reflect only the beam size and the signal-to-noise ratio of the 
SN detection and they are dwarfed by the systematic errors 
introduced by the atmosphere. We determine a best esti- 
mate for the radio SN position and associated uncertainty by 
calculating the weighted mean of the centroid positions and 
find a(J2000) = 02h43™27.9709s, (5(J2000) = +37°20'44.692" 
with uncertainties of 0.014 and 0.010 asec, respectively (la, 
standard deviation of the mean). 

The radio light curves are shown in Figure [T] spanning 
Af Ri 5 to 150 days. The observations reveal an early peak 
time of < 5 days at 8.46 GHz with an associated peak spec- 
tral luminosity of L^ w 9.7 x 10^^ erg s"' Hz"'. In compar- 
ison with the radio properties of other SNe Ibc, SN 2007gr 
is among the least luminous, about lO"' times lower than 
SN 19 98bw (associate d wifli GRB 980 425; Kul karni et al. I 
119981) and SN 200 9bb (ISod erberg et ^'2010) and most sim- 
ilar to SN 2002ap (iBerger. Kulkarni & Chevalier 2002). 

2.2. Chandra X-ray Observations 

We additionally observed SN 2007 gr with the CXO ACIS- 
I beg inning on 2007 A ug 29.1 UT (At « 16 days) for 20 
ksec (ISoderberg||2007h . We do not detect an X-ray source 
coincident with the optical and radio SN positions. Adopt- 
ing a power-law spectral model with photon index, F = 2, 
and a Galactic foreground column density of «h ~ 3.7 x 
10^" cm"^, we place a 3a upper limit on the X-ray flux of 
Fx < 3.9 X 10"'^ erg cm'^ s"' (0.2-10 keV). At the distance 
of the SN, this implies a luminosity of Lx ^ 4.0 x 10^^ erg s~' 
which is a factor of 10^ below the afterglow lu minosi- 
ties o f sub-energetic GRBs (e.g., GRB 980425; Plan et al. I 
I2OOOI) and similar to the X-ray luminosity of SN 2002ap on 
a comparable timescale dSutaria et al. 1 120()3|) . A compari- 
son with our nearly simultaneous VLA observations indi- 
cates a radio-to-X-ray spectral index steeper than sa -0.7, 
consistent with the measured values for other well-studied 
radio/X-ray SNe Ibc shortly after explosion and during the 
epoc h that inverse Compton em ission dominates the X-ray 
flux (IChevaUer & Franssonll2006l) . 

3. A SIMPLE MODEL FOR THE RADIO EMISSION 

The radio emission from SNe, including those associated 
with LGRBs, arises from the dynamical interact ion of the 
blast wave with the circumstellar medium (CSM; IChevalien 
1982; Sari. Piran & Naravanlll998h . In this process, flie blast- 
wave accelerates CSM electrons into a power-law distribu- 
tion, A^(7) oc 7"'\ above a minimum Lorentz factor, 7,,,. The 
interaction of the accelerated electrons with amplified mag- 
netic fields gives rise to non-thermal synchrotron emission. 
For SNe minimally affected by external absorption processes, 
a low frequency spectral turn-over is produced by synchrotron 
self-absorption (SS A) and defines the spectral peak frequency, 
i^p. In this scenario, the radio spectrum is given by /v oc i^^/^ 



below Vp and /v oc it'-P'^^I^ above Vp. As shown in Fig- 
ure |2] our multi-frequency radio observations of SN 2007gr 
are well-described by a synchrotron self-absorbed spectrum 
with /5 PS 3.2 across multiple epochs. We note that an extrapo- 
lation of the optically-thin synchrotron spectrum to the X-ray 
band is consistent with the observed CXO upper limit. 

3.1. Initial Constraints 

For radio SNe characterized by SSA, the radius of the blast- 
wave, R, and its time-averaged velocity, v, can be derived from 
the observed value of Vp and the peak sp ectral luminosity, L,y ^ 
(IChevahedll998l:lKulkarni et a/. 1119981) . For p « 3, the blast- 
wave radius is given by 
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where £« and €b are the relative fractions of post-shock en- 
ergy density shared by accelerate d electrons and amplifie d 
magnetic fields, respectively (Chevalier^ FranssonI |2006|). 
Throughout the paper we make the assumption that the radio 
emitting region is half of the total volume enclosed by a spher- 
ical blastwave. ForSN2007gr, wefindi-,, «4.0, 3.2, and 1.8 
GHz and U,,, « (9.8, 8.3, 8.3) x lO^^ erg s"' Hz"' at Ar « 
5, 8, and 18 days (Figure|2ll, and thus R « (4.0, 4.7, 8.3) x 
10'^ cm, respectively for 6^ = 6^. 

The minimum internal energy, £■„,!„, required to power 
the observed radio signal can be derived from the post- 
shock magnetic energy density, B^/^tt. As shown by 
IChevaher & FranssonI d2006.) . the amplified magnetic field is 
directly determined from the observed radio properties. 
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At Af w 5, 8, 18 days the SN 2007 gr radio spectra imply B ; 
0.55, 0.45, and 0.25 G. The minimum energy is 
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where, by definition of £■„,!„ we have assumed ee = ^b = 0.5. 
For our derived values of R and B across the three epochs, 
we infer Emin ~ (3.2, 3.5, 6.0) x 10"*^ erg. A more realistic 
scenario allows for a significant population of shocked pro- 
tons such that (e^ = e^) < 0.5. The total internal energy is then 
E = (0.5/eB) £™in and we find E w (1.6, 1.8, 3.0) x lO^*- erg 
fores = 0.1. 

In stark contrast, applying this same analysis to 
the observed radio emission from nearby GRB-SNe 
1998bw (Kulkarni e t al. 1998; Li & Chevalier 1999) and 
2006ai dSo derberg ej^ |2006|) as well as SN2009bb 
dSoderberg et al. . .20101) . points unambiguously to mildly- 
relativistic outflows with v ^ c and with energies of £ > 10"^^ 
erg. The properties of the SN2007gr blastwave are clearly 
dissimilar from those of relativistic explosions and are typical 
for ordinary, non-relativistic SNe Ibc. 

3.2. A Dynamical Model 

The temporal and spectral evolution of the radio emission 
can further constrain the properties of the SN 2007gr blast- 
wave and local environment. As the blastwave expands, the 



optical depth to SSA decreases and, in turn, Vp cascades to 
lower frequencies producing th e characteristic " bell-shaped" 
hght-curves of radio SNe Ibc (IChevalieiill998l). In model- 
ing thi s evolution, we adopt the formalism of lSoderberg et alA 
(2005) which assumes that the blastwave radius and the mag- 
netic field evolve in time as power laws R oc f "' and B oc t"" 
while £(, and es are const ant fractions of the post-shock energy 
density (IChevaliedll998l) . 

We perform a global fit of the multi-frequency light-curves 
for four parameters: Cf, Cr, ctr and as where Cf and Cr 
are normalization constants of the pea k flux density and opti- 
cal depth at a reference epoch, to (see ISoderberg et al. 1120051 
for a discussion of Cf and Cr). We find a reasonable fit 
for Cf « 3.7 x lO'^o g s°^ Cr « 3.0 x lO^^ s3■^ a, « 0.9 
and as ~ -1 at to = 10 days. The resulting model fits are 
shown in Figure[T]and imply /? « 5 x 10'^ {t/lO days)*^-^ cm 
^(i.e., V w 0.2c (f/10 days)-"') and B w 0.4 (f/10 days)"' G. 
Thus, the radio observations point to a freely-expanding SN 
in which the bulk ejecta travel with a constant velocity while 
the post-shock material is slightly decelerated. With these 
constraints, the the total internal energy of the radio emitting 
material is £ « 2 x 10'*'' (f/10 days)'" erg for e^ = eg = 0.1. 
These scalings are fully consistent with our preliminary esti- 
mates from individual-epoch analysis of the radio spectra in 

EH 

Next, the electron number density is n^ = (p- 2)/(p - 
1) (BV87r7„m,c2) w 2 x 10^ (r/5.5 x lO'^ cm)-^ cm'^ with 
7m « 2 (here we maintain the assumption that p « 3 and 
£(> = £« = 0.1). The density profile of the circumstellar en- 
vironment is thus consistent with the expectations for a stellar 
wind with a constant mass loss rate and wind speed. Assum- 
ing a typical Wolf-Rayet wind velocity o f v„ =10^ km s"' 
(e.g..'Capp a. Goss & van der Huchti[2004l) . and a nucleon-to- 
electron ratio of two (appropriate for a predominantly helium 
stellar wind), we infer a mass loss rate of M = 4TTni,mpR^v„ w 
6 X 10"^ M0 yr"' for the SN 2007 gr progenitor star. The tem- 
poral evolution of the physical parameters associated with the 
SN 2007gr blastwave and the circumstellar material are dis- 
played in Figure [3] We emphasize that the blastwave dynam- 
ics and physical properties that we deriv e for SN 2007gr are 
similar to those o f ordinary radio SNe Ibc (iBerger et al. 120031 : 
ISoderber^l2007l and references therein) an d fully consistent 
with dynamical models for radio SNe Ibc (IChevalieiJll998l : 
Chevalier & Fransson 2006). 

Finally we note that an independent constraint on the par- 
tition fractions can be obtained from the blastwave velocity 
since the observed synchrotron radiation requires a sufficient 
energy to accelerate electrons to relativistic speeds. For a non- 
or mildly-relativistic shock, the requirement that 7„, > 1 im- 
plies 
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which resul ts in a lower limit of e^ > 0.44 (v/O.lc)"^ fo r p ~ 3 
(Chevalier & Franssonl2006t see also lSoderberg et an i2005). 
For SN 2007gr, this constraint leads to e^ ^0.1 for v > 0.2c 
and is fully consistent with our equipartition model 03.2I) . 
therefore supporting our derived energy estimate. 

3.3. Bulk Ejecta Parameters 

B ased on mode l ing of the optical light-curves and spec- 
tra, iHunter et al. I (l2009l) reported values for the total ki- 
netic energy and mass of the SN ejecta of ^sn ~ (1-4) x 



lO''' erg and Mq « (2-3.5)M0, respectively. The veloc- 
ity of the bulk ejecta is thus Vbuik ~ (1-8 -Esn/^bj)'^^ ~ 
9,000-13,600 km s"' (Iwamoto et al. "20 03b, in line with 
the observed phot ospheric velocities ( Valenti et al. 1 120081; 
iHunter et al. 1120091) . Theoretical considerations predict that 
the coupling of energy and velocity within the homologous- 
expanding SN ejecta is characterized by 
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dMatzner & McKed 119991: iBerger. Kulkarni & Chevaheri 
I2OO2I) . For SN2007gr, the bulk ejecta parameters predict 
that the ejecta traveling at v > 0.2c carry an energy of 
E « (0.1-4.5) X lO^*" erg, which is fully consistent with 
the energy inferred from our free-expansion model ( ^3.2l l. 
Therefore there is sufficient energy in the high- velocity ejecta 
to account for the observed radio signal within the framework 
of a standard homologous explosion. 

4. A COMPARISON WITH PARAGIETAL. 2010 

As detailed in the previous section, the temporal and spec- 
tral evolution of the SN 2007 gr radio emission point to an or- 
dinary and non-relativistic blastwave with velocity, v w 0.2c, 
and energy, £ w 2 x 10"*^ erg, for e^ = e^ = 0.1. Similar blast- 
wave parameters are independently inferred from the bulk 
ejecta properties and exp ected energy profile of the ejecta 
(Eqn.|3l. In stark contrast. IParagi et al. I (120 lOl) propose a rel- 
ativistic blastwave velocity of v > 0.6c based on their analy- 
sis of a low signal-to-noise VLBI observation of the SN at 
Af K, 84 days from which they claim a lower limit on the 
blastwave radius, of 7? > 1 .3 x 10'^ cm. This is a factor of > 3 
larger than our modeling estimate. In the following sections 
we address this claim by first considering the implications of 
a mildly-relativistic outflow in the framework of our dynam- 
ical model for the SN 2007gr radio light-curves and spectra. 
Following this discussion, we present the results from our in- 
dependent analysis of the VLBI data. 

4.1. Severe Departures from Equipartition? 

As shown in Equation [T] the blastwave radius depends 

only weakly on the partition fractions, R oc (ee/ffi)''"^- 

This stems from the fact that the energy in electrons (£(,) 

and amplified magnetic fields {Eb) scale with the spherical 

1 1 o 

blastwave radius as Eg cc R and E^ o: R , respectively, 
such that the total energy budget (Eg + Eb) is minimized at 
equipartition (Figure |4]i. Therefore, in order to accommo- 
date a factor of > 3 increase in the time-averaged veloc- 
ity within the framework of our spherical blastwave model 
would require severe departures from equipartition, Ce/eB < 
10"^. Such deviations are unprecedented in astrophysical sys- 
tems and indeed, detailed VLBI studies of other radio SNe 
that point to relative partition fractions close to equiparti- 
tion ( e. g., £e/ea » 0.004 for SN 1993J; Fransson & Biornsson 
119981; IChandra. Ray & Bhatnagari 12004 ) while broadband 
modeUng of GRB afterglows typically indicates Cg/eB ~ 10 
(iPanaitescu & Kumadl2002l;lYost et aU Q.003). 

Furthermore, Eqn.|4] shows that a v > 0.6c outflow requires 
the relativistic electrons to harbor a significant fraction of 
the post-shock energy density, e^. > 0.01. This is a factor of 
10^ higher than the value required to accommodate the VLBI 
measurement within our spherical blastwave model. To rec- 



oncile this inconsistency would require an atypical modifica- 
tion of the electron energy distribution. 

Severe deviations from equipartition also impose a signif- 
icant increase in the total energy of the radio emitting ma- 
terial. In 33.21 we report the modest energetics required by 
a non-relativistic model for the blastwave in equipartition, 
£■ sa 2 X lO'*^ erg. However, a mildly-relativistic velocity of 
V > 0.6c would require a magnetically dominated blastwave 
with an enormous energy of E > 10^^ erg, exceeding the total 
bulk energy of the explosion (Figure IH. 



4.2. A Collimated Outflow? 

To a void severe departures from equipartition. lParagi et al. I 
(I2OIOI) propose that the outflow is highly collimated into jets 
with opening angles of just Oj w 15 degrees. Thus the ra- 
dio emitting region fills only a fraction of the total solid an- 
gle, fh = il-cosOj) « 0.03. Reducing the area of the ra- 
dio emitting region serves to increase the radius associated 
with equipartition, thus bringing the minimum of the (E^+Eb) 
curve closer to their claimed blastwave radius. In this sce- 
nario, they adopt a blastwave energy of £ « 3 x 10"*^ erg by 
allowing for modest departures from equipartition. The asso- 
ciated isotropic-equivalent energy would be, Ei^o = (E/ft) ~ 
10"*^ erg. iParagi et al. I (1201 Ol) further adopt a mass loss rate 
of M « 3 X 1 0~^ Mq , a factor of two lower than our own esti- 
mate. 

Such relativistic and highly collimated jets are are incom- 
patible with the s tandard expec tations for a homologous SN 
explosion. Thus, iParagi et al. I (1201 0) appeal to a model in 
which the blastwave detached from the bulk ejecta at the time 
of explosion, similar to the model for GRB-associated SNe. 
In this scenario, the trans-relativistic SN 2007 gr jets raced 
ahead of the bulk spherical outflow and freely-expanded^ until 
they swept up a circumstellar mass comparable to their own 
rest mass, caus ing a deceleration to non-relativis tic expansion 
at a time, fdec (IWaxman. Kulkarni & Fraillll998l) . Thereafter, 
the dynamics of the outflow approach the Sedov-Taylor so- 
lution characterized by a radial exp ansion, R oc f^/^, into a 
wind environment (e.g.. lSedovlll946h . Due to lateral spread- 
ing of the jets, the outflow also approaches spherical sym- 
metry within a dynamical timescale of fdec- This is roughly 
the time for the outflow radius to double in size, and thus 

tsph ~ ^ 'dec. 

Throughout its evolution leading up to fsph, the blast- 
wave experiences several dynamical phase transitions, from 
freely-expanding to decelerated expansion, relativistic to non- 
relativistic, and collimated to spherical evolution. These tran- 
sitions give rise to abrupt cha nges in the temporal decay 
of the radio light-curves (e.g.. Sari. Piran & HalpernI [19991 ; 



iFrail. Waxman & Kulkarnill2000l) . In particular, spherical and 
non-relativistic expansion requires a steepen ing of the opti- 
cally thin flux density to F^, oc f"^ ■' for /:> « 3 ('Chevali er & Lil 
I2OOOI) . In the case of SN 2007 gr, however, the optically thin 
flux evolution decays steadily throughout our radio observa- 
tions spanning At « 5- 150 days: F^ oc f-"-9±0-2 (L4 GHz), 
/V oc r''-^*"^ (4.9 GHz), and /v oc f-i"±oi (8.5 GHz; see Fig- 
ure[Tll. Thus we find no evidence for any phase transitions on 
this timescale, and we highlight that the observed temporal 
decay is statistically inconsistent with the expectations for a 
decelerated blastwave. 

^ Here we note that trans-relativistic outflows may not follow formal 
Blandford-McKee dynamics owing to their low bulk Lorentz factors. 



The timescale required for the jets to reach fdec depends 
on the rate at which the jets spread sideways. We con- 
sider two extreme cases that bracket the range of hydrody- 
namic evolutions: (i) lateral s preading with a rate, v ^ c, and 
(ii) m inimal spreading (see iGranot. Ramirez-Ruiz & Loebl 
l2005l for a full discussion). In the former case, /dec ~ 
100 (E/IO^' erg)(M/10-^^ Mq yr"')"' days, while the lat- 
ter case is longer by a factor of /^' since it assumes 
there is no lateral jet spreading durin g relativistic expan- 
sion (IChevaUer & Lill2000 : Waxman 2004). For Parag iefa/.l 
(l2010l) 's proposed jet parameters, a SN 2007 gr relativistic out- 
flow would decelerate at fdec ~ 1 - 30 days and approach 
spherical symmetry at fjph ~ 3-90 days. Therefore, even in 
the unlikely scenario that the SN 2007gr radio outflow was 
initially jetted, relativistic, and experienced minimal lateral 
spreading, by the epoch of the VLBI observation (Af w 84 
days) these ejecta would have roughly transitioned to non- 
relativistic and spherical expansion. Moreover, as noted 
above, the observed light-curves show no evidence for such 
phase transitions on these timescales. 

In summary, the exotic physical scenarios required by a 
mildly-r elativistic blastwav e with a velocity of v > 0.6c pro- 
posed bv'Paragi et al. I (12010) are overall inconsistent with the 
observed temporal and spectral evolution of the SN 2007gr ra- 
dio emission. We also show that a homologous free-expansion 
blastwave model with a non-relativistic velocity, v « 0.2c and 
an energy, £ « 2 x 10**^ erg, characterized by shock micro- 
physics near equipartition provide an excellent (and more nat- 
ural) description of the radio observations. We next turn to an 
independent analysis of the SN 2007gr VLBI dataset. 

4.3. Very Long Baseline Interferometry Data 

SN2007gr was observed with the European Very Long 
Baseline Interferometry Network (EVN) on two epochs be- 
ginning at 2007 Sep 6 and Nov 5 UT (Af « 24 and 84 days, 
respectively). We retrieved the data from the EVN archive'^ 
(programs RP007 and GP044; PI Pai-agi). Both EVN ob- 
servations were carried out at 4.9 GHz. The first observa- 
tion was carried out in real-time e-VLBI mode, lasted 12 
hours, and included the Darnhall, Medicina, Jodrell Bank, 
Onsala, Torun, and Westerbork (WSRT) antennas. The data 
were recorded with four 8 MHz sub-bands in dual polarization 
and 2-bit sampling, resulting in an aggregate data rate of 256 
Mbps. During the phase-referencing cycles, one minute was 
spent on the reference source J0253H-3835, and 4.5 minutes 
on SN 2007 gr After 2 cycles including SN 2007 gr, two ad- 
ditional background quasars (J0230H-4032 and J0247H-3254) 
were also observed while 3C 454.3 and 3C84 were used as 
fringe finders. 

The second observation included the Medicina, Jodrell 
Bank, WSRT, Cambridge, Torun, Noto, Onsala, Effelsberg, 
Hartebeestoek, and the Green Bank Telescope antennas. The 
observation lasted 10 hrs with a recording rate of 1 Gbps for 
the EVN stations (eight sub-bands of 16 GHz in dual polar- 
ization and 2-bit sampling), while the GBT observed with 512 
Mbps (same setup, but with 1-bit sampling). The observing 
scheme was similar to the first epoch, but the two background 
quasars were observed less frequently (only after 4 cycles on 
SN2007gr). 

The orig inal analysis of these data was reported by 
iParagi etal . (2010); here we present an independent analy- 
sis. The data reduction was performed using standard pack- 

' http://archive.jive.nl/scripts/portal.php 



ages within AIPS. Total electron content maps of the iono- 
sphere were used to correct for associated phase changes. A- 
priori amplitude calibration was applied using system temper- 
ature measurements and standard gain curves. We performed 
phase-calibration using the data for 3C 84 to remove instru- 
mental phase offsets among the frequency bands. We then 
fringe-fitted the data from J0253H-3835 and transferred the 
phase calibration to SN 2007gr 

4.3.1. VLBI Epoch 1 

In the first VLBI epoch, we confirm the detec- 
tion of an unresolved source at coordinates, a(J2000) = 
02H3™27.9715\ (5(J2000) = +37°20'44.687" (±5 mas in each 
coordinate, dominated by the positional uncertainty of the 
phase calibrator) which is within 2a of our weighted mean 
VLA SN position . We measure an integrated flux density 
of /v = 235 ± 60 /iJy for the source and note that this is the 
strongest source detected within the 3cr VLA localization re- 
gion (Figure |5]). The likelihood of detecting a source with 
S/N?s 3.9 (rms noise= 60 //Jy) in one one of the 26 indepen- 
dent VLBI beams within the 3(7 position ellipse is roughly 
3 X lO""'. This, along with the consistent flux of the source 
with our nearly coincident VLA measurement (/v.4,86 ghz = 
233 ±76 ^Jy at Sep 12.31 UT; Table [B supports the iden- 
tificati on of t he VL BI source as SN2007gr. We note that 
Paragi e7a/. I (1201 Ol) reported a higher integrated flux density 
for the SN by a factor of 1 .8 ± 0.6 while their associated rms 
noise level was also higher by a factor of 1.3. 

Since the SN emission is unresolved, the beam size (7.3 x 
6.7 mas) implies an upper limit on the diameter of the emitting 
region of < 9.3 x lO'^ cm. Assuming spherical expansion, 
this constrains the time-averaged velocity to be v < 7.4c. 

4.3.2. VLBI Epoch! 

In the second epoch, the longest baselines provided by 
the Green Bank Telescope and Hartebeesthoek antennas re- 
sult in a higher resolution image with a synthesized beam of 
3.0 X 0.93 ma s (natural weighting ). In their original analysis 
of these data, iParagi et al. I ( 1201 Oh reported the weak detec- 
tion of a source coincident with the first epoch SN position 
with flux density, /v ~ 60 /iJy (rms map noise, 13 /iJy). In 
comparison, coincident observations at j^obs = 4.9 GHz with 
the WSRT and VLA indicate some what higher SN flux den- 
sities of iv.4,9 GHz = 259 ± 40 ^Jy dParagi et al. I l20ly) ^ and 
146 ± 34 /iJy (Table [B, respectively 'Paragi et al. I (1201 0*) in- 
terpreted the WSRT- VLBI flux discrepancy to the radio emis- 
sion being resolved on the longest baselines. Adopting the 
high-resolution beam size as a lower limit on the angular di- 
ameter, they claimed that the apparently resolved SN emission 
(i.e., "missing flux") would require a time-averaged blastwave 
velocity of v > 0.6c. We note, however, that the VLA- VLBI 
flux discrepancy is less significant (consistent at the 2a level) 
and we further address the VLA- WSRT flux discrepancy at 
the end of this section in the context of host galaxy contami- 
nation. 

In our independent analysis of the second epoch VLBI data, 
we attempted to verify the source extension by excluding the 
data from the Green Bank and Hartebeesthoek antennas re- 
sulting in a lower resolution image with a synthesized beam 
size of 13 X 8.2 mas, comparable to that in the first VLBI 
epoch. We confirm the presence of a weak and apparently ex- 
tended source consistent with the position of the SN and mea- 
sure an integrated flux density of /v ~ 64 /iJy (rms noise level, 
(T = 13 /iJy per beam. Figure |6]l. Thus, a flux discrepancy 



persist s even on shorter baselines. We note that lParagief a/. I 
(1201 Ol) similarly reported no increase in the flux density when 
the data from the Hartebeesthoek and Greenbank antennas 
were excluded. If this effect is attributed to resolved emis- 
sion from a relativistic outflow, the required velocity is sig- 
nificantly higher, V > 2.6c, which implies far more stringent 
requi rements on the prop erties and dynamics of the blastwave 
than iParagi et al. I (l2010l) 's proposed value of v > 0.6c (see 

m. 

An alternative interpretation is that the observations suf- 
fer from a modest loss of phase coherence, a common 
phenomenon in phase-referenced VLBI observations which 
is difficult to correct at low signal-to-noise ratios. Such 
losses cause an apparen t decrease in the total flux density 
dMarti-Vidal et a/. 1 12010'). and could explain the modest flux 
discrepancy observed between the nearly coincident VLBI- 
VLA measurements"'. We therefore propose a more conser- 
vative interpretation in which the flux discrepancy and appar- 
ent extension of the SN - revealed even on shorter VLBI base- 
lines - are the result of modest phase decoherence affecting 
this low signal-to-noise observation. Under this interpreta- 
tion, we adopt the beam size as an upper limit on the apparent 
diameter of the otherwise unresolved SN emission to place 
a constraint on the blastwave velocity of v < 2.6c which is 
fully consistent with our non-relativistic spherical blastwave 
model. 

Finally, we consider the modest flux discrepancy observed 
for the VLA and WSRT measurements near the epoch of the 
second VLBI observation. Given the that the WSRT syn- 
thesized beam (3.2 asec) is significantly larger than that of 
our VLA observation (1.4 asec, B-array) we propose that the 
WSRT data may be contaminated by diffuse host galaxy emis- 
sion. To test this hypothesis, we analyzed VLA observations 
of the host galaxy, NGC 1058, obtained at 4.9 GHz prior to 
the discovery of SN2007gr and publicly available from the 
NRAO archive". We find that the SN is located in a strongly 
star-forming region of its host galaxy characterized by en- 
hanced radio emission. From these data, we measure diffuse 
host galaxy emission at the position of the SN with a flux 
density per beam of /v ~ 185 ± 15 /iJy (D-array, 15.4 asec 
beam) and /v « 23 ±29 /iJy (B-array, 1.4 asec beam). Given 
that the WSRT resolution is intermediate between these two 
cases, it is plausible that the WSRT observations are contam- 
inated by low-level and diffuse host emission. This is fur- 
ther supported by the fact that an earlier WSRT observation at 
Af w 24 days showed a similarly high flux density in compar- 
ison to our nearly contemporaneous VLA observations (see 
Figure [Til. 

In conclusion, we attribute the modest VLA- VLBI flux dis- 
crepancy and apparent extension of the SN emission in the 
low S/N second epoch VLBI observation to the likely loss of 
phase decoherence. We further attribute the VLA- WSRT flux 
discrepancy to contamination by diffuse host galaxy emission 
near the explosion site. 

5. CONCLUSIONS 

In this paper, we present a critical analysis of blastwave 
properties of Type I c SN 2007gr following the recent claim by 
iParagi et al. "(2010) that this otherwise ordinary SN produced 
a relativistic outflow similar to those of nearby LGRBs. We 
show that the full dataset of radio (VLA and VLBI) and X-ray 
observations for SN 2007 gr are more naturally explained by 



an ordinary, non-relativistic, and homologous SN explosion. 
Our conclusions stem from the fact that we directly measured 
the frequency and flux density of Vp across several epochs. 
This result underscores the necessity of multi-frequency and 
long-term radio monitoring of SNe Ibc in the search for rela- 
tivistic outflows. We conclude with the following points: 

1. A freely-expanding and non-relativistic blastwave 
model best reproduces our extensive radio and X-ray 
observations, indicating an expansion velocity of v w 
0.2c and energy of £ « 2 x lO"*^ erg. These blastwave 
parameters are consistent with those inferred from the 
bulk ejecta properties given the expected energy profile 
of the ejecta. 



2. iParagi et oTI ( l2010l) 's proposed mildly-relativistic 
blastwave velocity would require exotic physical 
scenarios (severe departures from equipartition or a 
decelerated, detached blastwave) that are implausible 
given the requirements of our VLA observations. 

3. Through our independent analysis of the VLBI data 
sets, we confirm the weak detection of a source con- 
sistent with our VLA position for SN 2007gr However, 
in our conservative interpretation of the low signal-to- 
noise second epoch detection, we attribute the apparent 
loss in flux density to modest phase decoherence, in- 
stead of relativistic SN expansion. We also suggest that 
the WSRT observation is contaminated by diffuse emis- 
sion from the host galaxy. 

Finally, we note that while the data for SN 2007gr do not 
point to a relativistic explosion, they do offer new insight 
on the non-thermal properties of Type Ibc SNe. In addi- 
tion to being one of the nearest of such explosions discov- 
ered to date, SN 2007 gr is also one of the least radio lumi- 
nous. This can be directly attributed to its low density cir- 
cumstellar environment that was shaped by a steady mass 
loss rate, M w 6 x 10"^ M© yr"' ( ^3.2l i. For comparison, 
this is 10"^ - 10^ times lower than the mass loss rates inferred 
for radio SNe Ibc with the st rongest circumstellar interac- 
tion (e.g., SNe 2003L, 2003bg; Soderberg et gLl|2Q0l ^200^) 
which are preferentially detected in the radio and X-ray bands 
thanks to their luminous non-thermal emission. Since SNe 
with low mass loss rates generally give rise to weak radio 
signals, they are only detectable nearby {d < 30 Mpc) with 
current cm-band facilities. This, together with the discovery 
rate of SNe Ibc within d k, \Q Mpc (a few each decade), sta- 
tistically suggests that SN 2007gr represents one of the most 
ordinary SNe Ibc studied to date, characterized by explosion 
and environmental properties that are typical of the bulk pop- 
ulation. With the significant improvement in continuu m sen- 
sitivit y enabled by the Expanded VLA (EVLA; Perlev e7a/. | 
l2009l) . we will soon be able to detect such ordinary SN Ibc 
to distances of ~ 100 Mpc which will, in turn, broaden our 
understanding of these unique cosmic explosions. 

The authors especially thank Mark Reid for helpful discus- 
sions. We also thank Dale Frail, Edo Berger, Andrew Mac- 
Fadyen and Eli Waxman. AMS is supported by a Hubble 
fellowship. RAC acknowledges support from NASA grant 
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'" Indeed. iParagi gf 011420101) similarly report that the data were "affected 
by modest phase eiTors". 

' ' http://archive.cv.nrao.edu/ 



REFERENCES 



Ai-nett, W. D. 1982, ApJ, 253, 785. 

Begelman, M. C. and Sarazin, C. L. 1986, ApJ, 302, L59. 

Berger, E., Kulkarni, S. R., and Chevaher, R. A. 2002, ApJ, 577, L5. 

Berger, E., Kulkarni, S. R., Frail, D. A., and Soderberg, A. M. 2003a, 

Astrophys. J., 599, 408. 
Berger, E. et al. 2003b, Nature, 426, 154. 

Cappa, C, Goss, W. M., and van der Hucht, K. A. 2004, AJ, 127, 2885. 
Chandra, R, Ray, A., and Bhatnagar, S. 2004, ApJ, 612, 974. 
Chevalier, R. A. 1982, Astrophys. J., 258, 790. 
Chevalier, R. A. 1998, ApJ, 499, 810. 
Chevalier, R. A. and Fransson, C. 2006, ApJ, 651, 381. 
Chevalier, R. A. and Li, Z. 2000, ApJ, 536, 195. 
Chomock, R., Filippenko, A. V., Li, W., Foley, R. J., Reitzel, D., and Rich, 

R. M. 2007, Central Bureau Electronic Telegrams, 1036, 1. 
Elias, J. H., Matthews, K., Neugebauer, G., and Persson, S. E. 1985, ApJ, 

296, 379. 
Filippenko, A. V. 1997, ARA&A, 35, 309. 
Filippenko, A. V. and Sargent, W. L. W. 1985, Nature, 316, 407. 
Frail, D. A., Soderberg, A. M., Kulkarni, S. R., Berger, E., Yost, S., Fox, 

D. W., and Hanison, F A. 2005, ApJ, 619, 994. 
Frail, D. A., Waxman, E., and Kulkarni, S. R. 2000, ApJ, 537, 191. 
Fransson, C. and Bjomsson, C. 1998, ApJ, 509, 861. 
Granot, J., Ramirez-Ruiz, E., and Loeb, A. 2005, ApJ, 618, 413. 
Hjorth, J. etal. 2003, Nature, 423, 847. 
Hunter, D. J. etal. 2009, A&A, 508, 371. 
Iwamoto, K., Nomoto, K., Mazzali, R A., and et al. 2003, in Supernovae 

and Gamma-Ray Bursters, ed. K. Weiler, volume 598 of Lecture Notes in 

Physics, Berlin Springer Verlag, 243. 
Kulkarni, S. R. et al. 1998, Nature, 395, 663. 
Li, Z. and Chevalier, R. A. 1999, ApJ, 526, 716. 
MacFadyen, A. 1., Woosley, S. E., and Heger, A. 2001, Astrophys. J., 550, 

410. 
Madison, D. and Li, W. 2007, Central Bureau Electronic Telegrams, 1034, 1. 
Marti-Vidal, I., Ros, E., Perez Torres, M. A., Guirado, J. C, Jimenez- 

Monfen-er, S., and Marcaide, J. M. 2010, ArXiv e-prints. 1003.2368 



Matzner, C. D. and McKee, C. F 1999, ApJ, 510, 379. 

Mazzali, P A. et al. 2003, ApJ, 599, L95. 

Panaitescu, A. and Kumar, R 2002, ApJ, 571, 779. 

Paragi, Z. et al. 2010, Nature, 463, 516. 

Perley, R. et al. 2009, IEEE Proceedings, 97, 1448. 

Plan, E. et al. 2000, ApJ, 536, 778. 

Podsiadlowski, P, Joss, P C, and Hsu, J. J. L. 1992, ApJ, 391, 246. 

Sari, R., Piran, T., and Halpem, J. P 1999, ApJ, 519, L17. 

Sari, R., Piran, T., and Narayan, R. 1998, ApJ, 497, L17-^. 

Sedov, L. I. 1946, Journal of Applied Mathematics and Mechanics, 10, 241. 

Silbermann, N. A. et al. 1996, ApJ, 470, 1. 

Soderberg, A. 2007a, The Astronomer's Telegram, 1205, 1. 

Soderberg, A. 2007b, The Astronomer's Telegram, 1187, 1. 

Soderberg, A. M. 2007c, Caltech PhD Thesis. 

Soderberg, A. M. et al. 2008, Nature, 453, 469. 

Soderberg, A. M. et al. 2010, Nature, 463, 513. 

Soderberg, A. M., Chevaher, R. A., Kulkarni, S. R., and Frail, D. A. 2006a, 

Astrophys. J., 651, 1005. 
Soderberg, A. M., Kulkarni, S. R., Berger, E., Chevaher, R. A., Frail, D. A., 

Fox, D. B., and Walker, R. C. 2005, Astrophys. J., 621, 908. 
Soderberg, A. M. et al. 2006b, Nature, 442, 1014. 
Soderberg, A. M., Nakar, E., Berger, E., and Kulkarni, S. R. 2006c, 

Astrophys. J., 638, 930. 
Sutaria, F K., Chandra, P, Bhatnagar, S., and Ray, A. 2003, A&A, 397, 1011. 
Tan, J. C, Matzner, C. D., and McKee, C. F 2001, Astrophys. J., 551, 946. 
Thompson, T. A., Chang, P, and Quataert, E. 2004, ApJ, 611, 380. 
Valenti, S. et al. 2008, ApJ, 673, L155. 
Waxman, E. 2004, ApJ, 602, 886. 

Waxman, E., Kulkarni, S. R., and Frail, D. A. 1998, ApJ, 497, 288. 
Wheeler, J. C. and Levreault, R. 1985, Astrophys. J.l, 294, L17. 
Woosley, S. E. and Bloom, J. S. 2006, Ann. Rev. Astr Ap., 44, 507. 
Yost, S. A., Hanison, F A., Sari, R., and Frail, D. A. 2003, ApJ, 597, 459. 









TABLE 1 
















VLA OBSERVATIONS OF SN 2007GR 






Date 


Fi,.lA 


FuA.9 


fi/.8.5 


F.,22 


Q2000 (8.5 GHz) 


52000 (8-5 GHz) 


Array 


(UT) 


(^Jy) 


(MJy) 


(/^Jy) 


(MJy) 


02''43"'+(s) 


+37''20'+(arsec) 


Config. 


2007 Aug 17.4 




777 ±52 


659 ±44 




27.9721 ±0.0007 


44.699±0.006 


A 


2007 Aug 18.4 


123 ±67 


954 ±60 


640 ±52 




27.9696±0.0008 


44.701±0.008 


A 


2007 Aug 19.5 








<327 






A 


2007 Aug 21.4 


142 ± 47 


687 ±58 


439 ±48 




27.9687±0.0011 


44.674±0.009 


A 


2007 Aug 24.5 


<228 


763 ± 64 


303 ± 56 




27.9718±0.0017 


44.653±0.018 


A 


2007 Aug 30.6 


677 ±42 


417 ±44 


172 ±42 




27.9751 ±0.0026 


44.696±0.021 


A 


2007 Sep 12.3 


628 ± 86 


233 ±76 


<218 












A 


2007 Sep 21.4 


478 ±159 


<238 














BnA 


2007 Sep 30.6 


529 ±77 
















BnA 


2007 Oct 23.4 


467 ± 60 


163 ±45 














B 


2007 Nov 18.3 


214 ±42 


146 ±34 














B 


2007 Dec 21.1 


< 146 


<157 














B 


2008 Jan 6.1 


145 ±27 
















B 



Note. — Upper limits are 3ct. 
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Fig. 1. — Our radio light-curves of SN2007gr as obtained with the VLA at Uabs = 1-4, 4.9, 8.5, and 22.5 GHz are well desc ribed by a freely-expanding 
blastwave model in which a non-relativistic shock accelerates electrons in a wind- stratified circumstellar environment (grey curves; j|3.2t . At 4.9 GHz, the VLB! 
flux densities at Af Ri 24 and 84 days are shown as encircled black dots (see j|4.3t . The WSRT flux density measurements at At ai 24 and 84 days (grey squares) 
are discrepant with our nearly simultaneous VLA measurements and could be due to contamination by underlying host galaxy emission. 
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Fig. 2. — The radio spectrum of SN2007gr across multiple epochs - A/ fs; 5.4 (red), 8.4 (green), and 17.6 (blue) days - is well described by a synchrotron 
self-absorbed spectral model with F oc u^l^ (Fy oc i/"<P"')/-) above (below) the spectral peak, Up. The observations indicate an electron energy index of /) fs 3.2. 
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Fig. 3. — The temporal and radial evolution of the physical parameters associated with the SN2007gr radio emitting material are shown in grey as derived 
from our dynamical model fit to the VLA data (see Figure[TJ. For comparison, we show the individual data points for the blastwave radius, energy, and magnetic 
field intensity as inferred from a single-epoch analysis of the SSA radio spectra at A? ss 5 (red), 8 (green) and 18 (blue) days (Figure|2). These individual points 
are overall consistent with our freely-expanding and non-relativistic dynamical model characterized by a wind-stratified circumstellar medium. 
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Fig. 4. — The total internal energy of the SN2007gr radio emitting material at At fs 84 days after explosion. Our modeling of the SSA radio spectrum 



across multiple ep ochs provides a tight constraint on the time-averaged velocity of the ejecta, v fs 0.2c, and a minimum internal energy of £,j 
at Af ; 



; 6 X 10"^ 



erg 



; 84 days (g3J) for e, 

147 



: eg = 0.5 (dashed black lines). Assuming these components each account for just 10 percent of the total internal energy, we find 
£ fa 2 X lO"" erg for e,, = eg = 0.1 (solid red line). Deviations from equipartition imply faster velocities for blastwaves energetically dominated by magnetic 
fields. Such deviations are also associated with a steep increase in the total energy budget. The mildly-relativistic blastwave properties reported bv lParagi et al. I 
(20101) (grey arrow) require extreme departures from equipartition, and assuming a spherical outflow implies a minimum internal energy of E si 10^^ erg for 
the blastwave (black encircled dot and an'ow). In this extreme scenario, the blastwave energy would rival the bulk explosion energy of £sN ~ (1 — 4) X lO" 
erg (horizontal grey bar) as estimated from modeling of the optical data (Hunter et al. 2009). As discussed in JI4.3I the VLBI observation suffers from a low 
signal-to-noise ratio and likely systematic effects which casts doubt on the claimed mildly relativistic expansion velocity. 
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Fig. 5. — EVN observation of SN2007gr on 2007 Sep 6 UT. The color scale extends from F^ = -250 to 270 /jjy per beam. The uncertainty associated with 
our weighted mean VLA position is shown as concentric error ellipses corresponding to la, 2a, and 3a. Within the la error ellipse, we detect an unresolved 
source with flux density Fi, = 235 it 60 fiJy which we identify as the SN (tick marks). 



14 



2007 Nov 6 UT 





w 

37 20 44.74 - 



44.72 - 



o 
o 
o 

CM 

^ 44.70 

O 



O 
lU 

a 



44.68 - 



44.66 m 



44.64 -- 



02 43 27.976 




27.974 



27.972 



27.970 



RIGHT ASCENSION (J2000) 



27.968 



Fig. 6. — EVN observation of SN 2007gr on 2007 Nov 5 UT, excluding the Green Banl; and Hailebeesthoek antennas. The color scale extends from 7v = —60 
to 60 fi]y per beam and the contours represent 1.5cr, 2a-, and 2.5cr levels (rms noise is 13/jJy per beam). At the position of the SN2007gr measured in the first 
VLBI epoch (black cross; la errors) we detect a weak and source with integrated flux density, F,, = 64 ± 3 1 ^Jy. Fitting a Gaussian in the image plane suggests 
that the emission region is apparently marginally extended with a major axis of 20+j(, mas, an upper limit of < 9.9 mas (3(t) for the minor axis, and a position 
angle of = 98;'; Jq degrees after convolution with the restoring beam. We attribute this apparent extension to the low signal-to-noise ratio of the image. 



